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Abstract 

Both  IrV  and  RhV  crystallize  in  the  adrV  structure,  with  a  transition  to  the 
higher  symmetry  Ll0  structure  at  high  temperature,  or  with  the  addition  of 
excess  Ir  or  Rh.  Here  we  present  evidence  that  this  transition  is  driven  by 
the  lowering  of  the  electronic  density  of  states  at  the  Fermi  level  of  the  cHrV 
structure.  The  transition  has  long  been  thought  to  be  second  order,  with 
a  simple  doubling  of  the  Ll0  unit  cell  due  to  an  unstable  phonon  at  the  R 
point  (0|f)-  We  use  first-principles  calculations  to  show  that  all  phonons  at 
the  R  point  are,  in  fact,  stable,  but  do  find  a  region  of  reciprocal  space  where 
the  Llo  structure  has  unstable  (imaginary  frequency)  phonons.  We  use  the 
frozen  phonon  method  to  examine  two  of  these  modes,  relaxing  the  structures 
associated  with  the  unstable  phonon  modes  to  obtain  new  structures  which 
are  lower  in  energy  than  Ll0  but  still  above  cHrV.  We  examine  the  phonon 
spectra  of  these  structures  as  well,  looking  for  instabilities,  and  find  further 
instabilities,  and  more  relaxed  structures,  all  of  which  have  energies  above 
the  cHrV  phase.  In  addition,  we  find  that  all  of  the  relaxed  structures,  stable 
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1.  Introduction 


A  major  goal  of  computational  condensed  matter  physics  is  the  determi¬ 
nation  of  structural  properties  for  compounds.  Once  the  structure  has  been 
determined,  other  properties,  e.g.,  strength,  ductility,  electronic  properties 
(including  superconductivity),  magnetism,  etc.,  can  be  determined.  Several 
methods  exist  to  attack  this  problem.  One  possibility  is  to  exhaustively 
search  an  experimental  database  of  known  structures  [it],  determining  the  low 
energy  structures  for  each  composition  of  the  target  materials.  Other  meth¬ 
ods  use  first-principles  calculations  on  a  small  set  of  target  structures  to  de¬ 
termine  parameters  which  can  be  used  to  predict  properties  of  more  complex 
systems.  These  include  the  Cluster  Expansion  Method j3,  0,  3,  di  0], tight- 
binding  paranretrization  methods [3-  §j,  and  atomistic  potential  methods^. 
Methods  of  the  first  type  can  be  combined  with  methods  of  the  second  type 
for  more  through  searches [lOi] . 

These  programs,  can,  however,  only  be  implemented  with  a  thorough 
knowledge  of  the  behavior  of  compounds.  Exhaustive  searches  of  a  database 
require  an  extensive  database  to  search,  including  all  structures  which  are 
known,  or  thought  likely,  to  form  for  a  target  system.  Parameterized  methods 
may  not  be  able  to  reach  all  regions  of  phase  space,  and  also  need  a  database 
of  structures  for  testingfllL  12] .  For  both  methods,  then,  it  is  useful  to  look 
at  less-common  structures  found  in  nature. 

The  adrV  structure,  shown  in  Figure  HJ  is  found  in  nature  only  in  the  pro¬ 
totype  compound |l3l.  114]  and  its  neighbor  in  the  periodic  table,  RhV  15l.  Il6  . 
This  orthorhombic  structure,  space  group  Cmmm  (#65),  can  be  viewed  as 
a  doubled  unit  cell  distortion  of  the  CsCl  structure [171]  or  the  tetragonal  Ll0 
structure,  space  group  P4/mmm  (#123)  [r|.  Indeed,  both  the  IrV  and  RhV 
phase  diagrams [19|  show  the  Ll0  structure  as  the  ground  state  for  Vanadium- 
poor  (x  =  40  —  48%)  [Ir,Rh]1_3,Va;,  and  Vanadium- rich  adrV  is  known  to 
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Figure  1:  The  alrV  structure.  The  box  indicates  the  boundaries  of  the  full  orthorhombic 
unit  cell.  In  the  Llo  or  (CsCl)  structure  the  Iridium  atoms  would  be  in  a  tetragonal  (cubic) 
arrangement  around  each  of  the  Vanadium  atoms. 

transform  to  Ll0  at  temperatures  above  506°  cfl8|. 

The  seemingly  straightforward  transition  pathway  from  Ll0  to  cJrV  was 
analyzed  by  Chen  and  Franzen(l8l|  in  the  context  of  Landau  theory:  First, 
double  the  tetragonal  Ll0  unit  cell  along  the  y  and  z  directions,  corresponding 
to  a  phonon  at  the  R  (0||)  point  in  reciprocal  space.  Second,  construct  the 
primitive  base-centered  orthogonal  unit  cell,  and  allow  the  two  Ir  atoms  in 
the  cell  to  move  opposite  each  other  along  the  z-axis,  while  the  V  atoms 
move  along  the  y- axis.  The  resulting  structure  is  symmetrically-equivalent 
to  that  of  adrV  and,  assuming  a  second-order  phase  transition,  relaxes  into 
the  adrV  structure. 

This  simple  picture  is  not  quite  accurate.  A  true  second  order  transi¬ 
tion,  as  described  by  Chen  and  Franzen,  requires  a  continuous  energy  path 
from  Ll0  to  adrV.  Thus  the  energy  of  the  Ll0  structure  would  be  lowered 
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by  making  an  infinitesimal  displacement  of  the  type  described  above,  cor¬ 
responding  to  a  phonon  instability  at  the  R  point.  We  have  performed 
first-principles  density  functional  theory  calculations  for  the  phonon  frequen¬ 
cies  at  the  R  point,  both  by  the  frozen  phonon  method  0  and  by  linear 
23l| .  We  find  that  all  of  the  modes  here  have  real  frequen- 


response[21l.  1221 
cies,  i.e. ,  they  do  not  lead  to  an  instability. 

Thus  the  transition  from  the  Llo  phase  to  the  adrV  phase  does  not  pro¬ 
ceed  through  the  simple  unit-cell  doubling  picture  described  above.  Instead, 
as  we  shall  show,  the  Ll0  phase  has  vibrational  instabilities  in  another  re¬ 
gion  of  the  Brillouin  zone.  This  paper  will  discuss  the  structures  arising  from 
these  instabilities,  and  explore  the  possibility  that  a  second-order  transition 
might  go  through  one  of  those  phases. 

In  Section  [2]  we  show  the  results  of  density  functional  theory  calculations 
for  the  energy  and  electronic  structure  of  the  cdrV  and  Ll0  phases  of  IrV, 
RhV,  and  some  neighboring  compounds.  We  show  that  in  IrV  and  RhV  the 
Llo  structure  has  a  relatively  high  electronic  density  of  states  at  the  Fermi 
level,  and  so  the  transition  to  cdrV  has  Jahn- Teller  character [24|. 

I11  Section  [3]  we  look  at  the  phonon  spectra  of  the  cdrV  and  Llo  phase. 
Not  surprisingly,  we  find  that  the  cdrV  phase  has  no  phonon  instabilities. 
We  do,  however,  find  that  the  Ll0  structure  is  vibrationally  unstable  on  and 
near  the  line  (a^|)  in  reciprocal  space. 

I11  Section  [4]  we  use  the  frozen-phonon  method  on  the  unstable  phonon 
modes  to  search  for  new  structures  which  have  lower  energy  than  Llo-  Such 
a  search  is  computationally  bound,  so  we  looked  at  all  unit  cells  with  eight 
atoms  or  less,  and  two  structures  with  32  atom  unit  cells.  In  this  range  we  find 
no  instability  which  relaxes  to  the  a  IrV  phase,  however  it  is  not  impossible 
that  searching  through  larger  unit  cells  would  find  such  an  instability,  in 
which  case  the  transition  would  indeed  be  second  order.  We  do  find  several 
new  structures,  some  of  which  have  apparently  never  been  seen  in  nature, 


and  one  which  was  previously  known  16 


In  Section  [5]  we  discuss  our  results,  and  thoughts  on  the  order  of  the 
Llo- cdrV  transition. 


2.  Energetics  and  Electronic  Structure  of  the  Ll0  and  alrV  Phases 

All  computations  were  made  using  the  Kohn-Sham{25[  formulation  of 
Density  Functional  Theory [26|  with  the  Perdew-Burke-Ernzerhofj27|  gener¬ 
alized  gradient  approximation.  Depending  on  our  needs,  we  used  the  Vienna 
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Figure  2:  Energy /volume  curves  for  IrV  and  RhV  in  the  Llo  and  olrV  structures,  deter¬ 
mined  from  VASP  and  QE  calculations  as  discussed  in  the  text.  For  ease  of  comparison 
we  have  set  the  minimum  energy  of  the  LIq  phase  to  zero. 


ab  initio  Simulation  Package  (VASP) 


30]  with  projector  augmented- 


wave  (PAW)  potentials [3 II] .  or  the  Quantum  Espresso  (QE)  package [321]  with 
the  supplied  ultra-soft  pseudopotentials  for  Ir,  Rh,  and  V. 

We  used  rather  large  plane  wave  cutoffs  of  350  eV  in  VASP  and  540  eV  in 
QE  to  ensure  convergence.  Similarly,  we  used  12  x  12  x  8  T-centered  k-point 
meshes  in  both  the  Ll0  and  adrV  structures,  leading  to  140  and  215  points, 
respectively,  in  the  irreducible  Brillouin  zones.  We  summed  over  electronic 
states  using  a  Fermi-Dirac  distribution [33j  with  a  temperature  of  65  meV 
(0.005  Rydbergs).  These  values  give  well- converged  results. 

We  show  our  results  for  IrV  and  RhV  in  Table  |T]  and  Fig.  [2]  In  both 
cases  we  see  that  the  VASP  PAW  potentials  and  the  QE  ultrasoft  pseudopo¬ 
tentials  are  in  excellent  agreement  with  one  another,  and  in  good  agreement 
with  experiment,  within  the  usual  errors  of  Density  Functional  theory  in  the 
generalized  gradient  approximation.  In  both  compounds  the  adrV  state  is 
below  the  Ll0  state  by  approximately  55  meV/formula  unit,  in  agreement 
with  experiment. 

From  experiment [isl]  we  know  that  Ll0  is  the  preferred  high  temperature 
structure  for  olrV.  Using  the  COMSUBS  routine  from  the  ISOTROPY [dl. 


35]  package,  we  find  that  an  orthorhombic  distortion  of  the  LIq  structure 


reduces  the  symmetry  from  space  group  P4/mmm  to  Cmmrn,  the  space 
group  of  the  adrV  structure.  It  is  plausible  to  argue  [r|  that  this  lower¬ 
ing  of  symmetry  is  the  pathway  for  the  Ll0  — >  alrV  phase  transition.  The 
symmetry-breaking  character  of  the  transition  is  evident  from  Fig.  [3l  where 
we  plot  the  electronic  density  of  states  of  both  phases  near  the  Fermi  level. 
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Table  1:  Equilibrium  lattice  constants  (in  A),  atomic  positions,  and  equilibrium  bulk 
modulus  (K0,  in  GPa)  for  the  Ll0  and  alrV  structures  of  IrV  and  RhV  determined  from 
experiment [131  Il5j| .  VASP  PAW  calculations,  and  QE  ultrasoft  pseudopotential  calcula¬ 
tions.  Note  that  we  give  the  primitive  tetragonal  lattice  parameters  for  the  Llo  structure, 
rather  than  the  more-common  face-centered  tetragonal  setting,  so  that  when  c/a  =  1  the 
Ll0  structure  reduces  to  the  cubic  CsCl  structure.  The  final  row  of  the  table  shows  the 
energy  difference  between  the  two  phases  (meV/formula  unit). 


IrV _ 

Ll0  Structure 


RhV 


Exp. 

VASP 

QE 

Exp. 

VASP 

QE 

a 

2.749 

2.755 

2.762 

2.754 

2.739 

2.745 

c 

3.651 

3.677 

3.668 

3.599 

3.660 

3.648 

K0 

272 

277 

230 

ct  IrV  Structure 

a 

5.791 

5.838 

5.816 

5.78 

5.849 

5.813 

b 

6.756 

6.759 

6.767 

6.65 

6.707 

6.725 

c 

2.796 

2.814 

2.823 

2.78 

2.792 

2.802 

K0 

271 

276 

227 

Ir/Rh  (4j) 

0.22 

0.216 

0.216 

0.214 

0.215 

V  (4g) 

0.28 

0.296 

0.297 

0.296 

0.297 

A  E 

53.6 

67.0 

50.4 

59.7 
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Figure  3:  Electronic  density  of  states  for  the  Ll0  (red)  and  cdrV  (blue)  phases  of  IrV 
(left)  and  RhV  (right),  found  by  smearing  the  eigenvalues  computed  by  VASP  using  a 
Fermi-Dirac  distribution  at  T  =  65  meV  (5  mRYl .  [331] 


As  we  are  only  interested  in  the  overall  behavior  of  the  density  of  states,  we 
compute  these  curves  by  smearing  out  each  eigenvalue  found  by  VASP  using 
a  Fermi  distribution  at  a  temperature  of  5  mRy.  We  see  that  at  the  Fermi 
level  the  density  of  states  is  twice  as  large  in  the  Ll0  phase  as  it  is  in  the 
adrV  phase.  This  is  pxmsistent  with  a  Jahn- Teller-like  symmetry  breaking 
and  phase  transition  24].  As  an  aside,  we  also  note  that  there  is  a  minimum 


in  the  Ll0  density  of  states  just  above  the  Fermi  level.  This  is  consistent  with 


the  phase  diagram  of  IrV  19],  which  shows  that  Iridium-rich  IrV  has  the  Llc 


structure.  Assuming  the  additional  Ir  replaces  V  in  the  Ll0  structure,  and 
using  a  rigid-band  model,  we  can  see  how  adding  Ir  to  the  system  would  raise 
the  Fermi  level  and  lower  the  density  of  states,  leading  to  a  more  stable  Llo 
phase.  We  find  a  similar,  though  less  pronounced,  state  of  affairs  in  RhV,  as 
is  also  seen  in  Fig.  |3] 

We  investigated  the  possibility  that  the  adrV  structure  might  be  found 
in  compounds  neighboring  IrV  and  RhV  in  the  periodic  table.  The  results 
are  shown  in  Table  [2]  Except  for  CoPt,  all  of  the  compounds  shown  here 
have  14  electrons  in  the  valence  band,  just  as  in  IrV  and  RhV.  In  all  cases 
except  CoPt  we  found  that  the  adrV  structure  was  locally  stable,  that  is, 
if  we  started  from  the  lattice  parameters  of  the  adrV,  replacing  Ir  and  V 
by  the  indicated  atoms,  and  allowed  the  system  to  relax,  it  remained  in  the 
adrV  structure  and  did  not  relax  to  the  higher  symmetry  Llo  structure.  This 
suggests  that  the  formation  of  the  adrV  structure  is  a  result  of  a  Fermi  level 
effect  associated  with  14  electrons  in  a  tetragonal  unit  cell.  We  did  not  check 
the  elastic  and  vibrational  stability  of  the  resulting  structure.  However  the 


7 


Table  2:  Energy  difference,  in  meV  per  formula  unit  between  the  Llo  and  odrV  structures 
for  the  given  compounds,  as  determined  by  VASP/PAW/PBE  calculations.  A  positive 
number  indicates  that  the  Llo  structure  is  lower  in  energy.  An  energy  difference  of  zero 
indicates  that  the  odrV  structure  relaxes  into  the  Llo  structure.  Ground  state  structures 
are  from  the  ASM  Online  Handbook[  19|. 

Compound  CoPt  CoV  IrNb  IrTa  IrV  RliNb  RliTa  RhY 

Ground  State  Ll0  D8&  Ll0  ctlrV  Ll0  odrV 

AE  0  2573  597o  79.9  -53.6  50T  706  -50.4 


Figure  4:  Electronic  density  of  states  of  the  ground-state  Ll0  and  hypothetical  cdrV  phases 
of  IrNb,  found  by  smearing  the  eigenvalues  computed  by  VASP. 


Llo  structure  was  lower  in  energy  for  all  compounds  except  IrV  and  RhV. 
Also,  as  shown  in  Fig.  [4l  for  these  compounds  the  density  of  states  of  the  Ll0 
structure  is  very  near  to  that  of  the  ctlrV  structure,  precluding  a  Jahn- Teller 
lowering  of  the  energy  with  decreased  symmetry.  We  therefore  conclude  that 
it  is  unlikely  to  find  the  odrV  structure  anywhere  except  in  IrV  and  RhV. 


3.  Stability  of  the  Llo  and  alrV  Phases 

A  system  is  stable,  or  metastable,  only  if  it  has  real  phonon  frequen¬ 
cies  for  all  k-points  and  modes  except  the  acoustic  phonons  at  the  T  point, 
which  are  guaranteed  to  vanish  because  of  translational  symmetry.  Since  the 
elastic  constants  are  related  to  the  long-wavelength  behavior  of  the  acoustic 
phonons[37}  this  implicitly  includes  the  Born  stability  criteria|38i] .  This  is  for- 


Figure  5:  The  base-centered  orthorhombic  alrV  structure.  Left:  Phonon  frequencies  IrV  at 
equilibrium,  determined  using  the  linear  response  method  with  ultrasoft  pseudopotentials 
from  the  Quantum  Espresso  package.  Right:  First  Brillouin  zone.  Symmetry  lines  are 
described  by  the  notation  of  Miller  and  Love  36j. 


tunate,  since  the  orthorhombic  ctlrV  structure  would  require  us  to  compute 
nine  elastic  constants. 

To  check  the  stability  of  the  adrV  phase  and  the  possible  meta-stability 
of  the  Ll0  phase  we  computed  phonon  frequencies  throughout  the  respective 
Brillouin  zones  using  linear  responsejill.  22l.  23|.  as  implemented  in  the  Quan¬ 
tum  Espresso  package  PHon  code|32|] .  Phonon  frequencies  were  computed  on 
a  reciprocal  space  grid  (“q-points”)  using  an  8  x  8  x  4  mesh  for  the  Llo  struc¬ 
ture  (45  points  in  the  irreducible  Brillouin  zone)  and  a  6  x  6  x  4  mesh  for 
the  a  IrV  structure  (39  points).  Linear  response  calculations  at  these  points 
yield  a  reciprocal-space  dynamical  matrix.  These  matrices  are  transformed 
into  a  real-space  dynamical  matrix  which  can  be  used  to  compute  phonon 
frequencies  over  the  entire  Brillouin  zone. 

This  approach  can  lead  to  aliasing  at  points  off  of  the  q-point  mesh.  To 
check  that  this  does  not  occur  we  also  performed  frozen-phonon  calculations (2oj] . 
wherein  we  construct  a  superccll  commensurate  with  the  given  q-point  and 
measure  the  change  in  energy  as  a  function  of  atomic  displacement  within  the 
supercell.  To  do  this  we  used  the  program  FROZSL,  part  of  the  ISOTROPY 


package [34 135] .  with  atomic  displacements  of  0.1  Bohr.  Electronic  structure 


calculations  were  then  performed  with  the  same  energy  cutoff  as  the  original 
unit  cell,  and  when  possible  the  same  k-point  mesh,  albeit  folded  back  into 
the  smaller  Brillouin  zone  associated  with  the  supercell. 

In  the  following,  in  interest  of  saving  space,  we  only  show  the  results  for 
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Figure  6:  The  tetragonal  Llo  structure.  Left:  Phonon  frequencies  for  IrV  at  its  minimum 
energy  structure,  determined  using  the  linear  response  method,  using  ultrasoft  pseudopo¬ 
tentials  from  the  Quantum  Espresso  package.  Frequencies  below  zero  (and  in  red)  are 
unstable  modes  which  actually  have  imaginary  frequencies.  The  dots  represent  calcula¬ 
tions  using  the  frozen  phonon  method  with  displacements  computed  by  FROZSL.  The 
second  panel  shows  the  frequencies  along  the  line  ( ).  Right:  First  Brillouin  zone,  with 
symmetry  lines  labeled  according  to  the  notation  of  Miller  and  Love^dj.  The  arrow  from 
U  to  S  indicates  the  region  of  unstable  phonons. 


IrV.  Our  calculations  for  RhV  show  a  similar  pattern. 

Fig.  [5]  shows  the  phonon  frequencies  of  the  equilibrium  structure  of  adrV. 
High  symmetry  points  and  lines  are  labeled  according  to  the  convention  of 
Miller  and  Love  3^].  All  the  phonons  have  real  frequencies,  confirming  that 
adrV  is  at  least  a  metastable  structure  for  IrV,  as  suspected.  There  is  no 
evidence  of  any  aliasing  in  these  calculations. 

Figure  [6]  shows  the  phonon  frequencies  of  the  Ll0  phase  of  IrV,  along 
symmetry  lines  labeled  according  to  the  convention  of  Miller  and  Love  [36| . 
We  show  results  of  frozen  phonon  calculations,  showing  good  agreement  be¬ 
tween  the  two  techniques,  even  off  the  linear  response  q-mesh.  As  expected, 
there  are  regions  of  reciprocal  space  with  imaginary  phonon  frequencies.  Sur¬ 


prisingly,  this  region  is  not  near  the  R  point  (O^)  as  described  by  Chen  and 
Fr arizen [l8|,  but  along  and  near  a  line  from  the  midpoint  of  the  U  line  to  the 
midpoint  of  the  S  line,  as  shown  in  more  detail  in  the  second  panel  of  the 
figure.  We  will  discuss  the  implications  of  this  instability  in  the  next  section. 


4.  Searching  for  Low  Energy  Structures 

The  previous  section  showed  that  the  adrV  structure  is  the  ground  state 
(or,  at  least,  a  low-energy  metastable  state)  of  IrV,  while  the  Ll0  structure  is 
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unstable  to  phonons  along  and  near  the  line  in  the  Brillouin  zone,  but 

not  at  the  R  point,  where  the  phonon  modes  are  real,  invalidating  the  second- 
order  phase  transition  scenario  of  Chen  and  Fr arizen [lsj] .  In  this  section  we 
will  look  to  see  if  it  is  possible  to  find  another  continuous  transition  path 
from  Ll0  to  «IrV. 

In  the  frozen  phonon  method,  we  displace  atoms  from  their  equilibrium 
sites  in  a  superccll  consistent  with  the  phonon  wave  vector  and  in  directions 
which  maintain  the  symmetry  of  the  phonon  mode.  The  frequency  of  the 
mode  is  then  directly  related  to  the  square  root  of  the  curvature  of  the  energy 
as  a  function  of  atom  displacement.  A  mode  with  an  imaginary  frequency 
will  then  as  a  matter  of  course  be  related  to  a  superccll  calculation  which 
has  a  negative  curvature,  leading  to  supercells  with  energy  lower  than  the 
original  state.  These  structures  can  then  be  relaxed,  leading  to  new,  or  at 
least  different,  structures. 

As  an  example  of  this,  consider  the  unstable  Ll0  phonon  at  the  point 
(\\\)  on  the  high  symmetry  “S”  line.  The  linear  response  calculations  find 
this  phonon  frequency  to  be  81*  cm-1.  To  generate  a  superccll  for  this  cal¬ 
culation  we  used  the  FROZSL  code  from  the  ISOTROPY  package [34.!.  f35j] . 
This  package  tells  us  that  this  point  has  phonons  in  three  irreducible  repre¬ 
sentations,  each  with  two  associated  modes.  For  each  representation  we  run 
three  total  energy  calculations  for  displacements  of  the  atoms  in  different 
directions,  with  the  maximum  displacement  of  0.1  Bohr.  The  energy  differ¬ 
ences  between  these  structures  and  the  ground  state  are  used  to  determine 
the  dynamical  matrix  at  this  point,  and  the  mode  frequencies. 

For  the  equilibrium  Ll0  parameters  found  in  Table  [D  the  S3  irreducible 
representation  (in  the  notation  of  Miller  and  Love[36|)  has  one  mode  with  an 
imaginary  frequency  of  74  cm-1.  (The  discrepancy  between  linear  response 
and  frozen  phonon  calculations  is  due  to  the  use  of  different  k-point  meshes 
and  the  anharmonicity  of  the  mode.  For  our  purposes  we  are  only  interested 
in  showing  that  the  mode  is  unstable  in  both  cases,  and  so  will  not  try  to 
refine  the  calculations  to  improve  the  agreement  between  the  two.) 

Diagonalizing  the  dynamical  matrix  allows  us  to  find  a  supercell  with 
displaced  atoms  which  has  an  energy  lower  than  the  Llo  phase.  This  superccll 
has  space  group  Cmmm  (^65),  and  with  the  appropriate  choice  of  origin 
has  Ir  atoms  at  the  (4e)  and  (4h)  Wyckoff  positions,  and  V  atoms  at  the 
(2a),  (2b),  and  (4j)  positions.  This  is  crystallographically  equivalent  to  the 
intermetallic  GasPts  structure  [39j,  and  we  will  refer  to  it  as  such.  Upon 
relaxation  we  find  a  minimum  energy  structure  which  has  approximately  the 
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Figure  7:  Linear-response  phonon  spectra  of  two  eight-atom  unit  cells  derived  from  un¬ 
stable  phonon  modes  of  Llo  IrV.  Left:  The  metastable  RI15V3  structure.  Right:  The 
Amm2  structure  with  Ir  atoms  at  two  (4d)  Wyckoff  positions  and  V  atoms  at  four  (2b) 
Wyckoff  positions.  The  imaginary  frequencies  of  the  unstable  modes  are  shown  as  negative 
frequencies  and  highlighted  in  red. 


same  density  as,  and  an  energy  20.5  meV/formula  unit  below,  the  relaxed 
Llo  structure  of  Table  [1]  Note  that  this  is  still  well  above  the  energy  of  the 
adrV  structure. 

In  a  similar  fashion,  if  we  look  at  the  phonons  at  0",  on  the  “U”  line,  we 
find  an  imaginary  U3  mode  with  a  frequency  of  70  cm-1.  This  mode  leads  to 
another  superccll  with  space  group  Cmnun,  but  now  the  Ir  atoms  are  on  the 
(4e)  and  (4g)  Wyckoff  sites,  while  the  V  atoms  occupy  a  pair  of  (4j)  Wyckoff 
sites.  The  relaxed  structure  again  has  a  density  comparable  to  Ll0,  but  its 
energy  is  only  11.0  meV/formula  unit  below  the  Ll0  structure  minimum.  In 
the  discussion  below  will  refer  to  this  structure  as  Cmnun. 

Neither  the  GaaPts  nor  the  Cmnun  structure  will  relax  to  the  ground 
state  adrV  structure.  However,  both  structures  have  imaginary  frequency 
long-wavelength  optical  mode  phonons.  If  we  venture  away  from  the  T  point 
we  will  have  to  deal  with  frozen-phonon  calculations  on  unit  cells  having 
hundreds  of  atoms,  so  we  will  only  examine  the  unstable  modes  at  T  in  both 
structures. 

First  consider  the  imaginary  frequency  Tj  mode  of  the  GasPts  struc¬ 
ture.  It  is  associated  with  a  superccll  of  space  group  Amm2  (#38),  with 
Ir  atoms  at  the  (4d)  and  (4e)  Wyckoff  sites  and  V  atoms  on  pairs  of  (2a) 
and  (2b)  sites.  This  structure  is  crystallographically  equivalent  to  the  Rh5V3 
structure  fit!  40']  and  will  be  referred  to  by  that  name.  The  relaxed  struc¬ 
ture  is  41.1  meV/formula  unit  below  the  Llo  structure,  with  similar  density. 
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This  structure  is  metastable,  as  seen  in  the  linear-response  phonon  spectrum 
plotted  on  the  left-hand  side  of  Fig.  0 

The  second  structure,  derived  from  the  unstable  TiT  mode  of  Cmmm, 
again  with  space  group  Amm2  (#38),  puts  the  Ir  atoms  on  two  (4d)  sites 
and  the  V  atoms  on  four  distinct  (2b)  sites.  The  relaxed  structure  has  an 
energy  24.1  meV/formula  unit  below  the  Llo  structure,  again  with  a  similar 
density.  This  structure  has  not  been  found  in  the  intermetallic  literature.  It 
is  vibrationally  unstable,  as  can  be  seen  from  the  phonon  spectrum  on  the 
right-hand  side  of  Fig.  0 

The  four  structures  Ga3Pts,  Rh5V3,  Cmmm  and  Amm2  are  the  only 
structures  with  eight  or  fewer  atoms  in  the  primitive  cell  which  can  be  gener¬ 
ated  by  relaxing  the  unstable  phonon  modes  of  the  Llo  structure.  If  we  wish 
to  examine  larger  unit  cells,  we  can  continue  along  this  line  of  research  ad  in¬ 
finitum,  or  until  we  find  a  supercell  where  the  atoms  relax  to  the  ground  state 
cdrV  structure,  a  supercell  with  energy  below  the  ground  state  structure,  or 
we  run  out  of  computational  resources.  In  most  cases,  including  this  one,  we 
will  reach  the  latter  limit  first.  We  did  compute  the  relaxed  energies  of  two 
supercells  associated  with  unstable  phonon  modes  in  the  Amm2  structure, 
and  one  with  an  unstable  mode  of  the  Cmmm  structure.  These  structures 
had  thirty-two  atoms  in  the  supercell,  and  energies  below  the  parent  struc¬ 
ture  but  above  the  metastable  Rh5V3  structure. 

Tables  0  0  and 0  summarize  all  of  the  calculations  discussed  here,  includ¬ 
ing  the  structural  derivation,  space  group,  lattice  constants,  atomic  positions, 
energy,  and  stability. 

5.  Discussion 

We  have  shown  by  direct  calculation  that  the  transition  from  the  Llo  to 
the  olrV  structure  of  both  IrV  and  RhV  is  a  result  of  a  Jahn- Teller  driven 
distortion  of  the  high-symmetry  unit  cell.  This  distortion  does  not  result 
from  a  zone-doubling  unstable  phonon  at  the  R  point  of  the  Ll0  Brillouin 
zone,  and  so  the  simple  Landau  theory  picture  does  not  hold.  This  does  not 
completely  eliminate  the  possibility  that  the  transition  is  second  order,  for 
we  found  an  entire  region  of  reciprocal  space,  in  and  around  the  line  (x||), 
where  the  Llo  structure  has  vibrational  instabilities.  We  examined  several 
of  those  instabilities,  and  found  that  they  lead  to  numerous  new  structures, 
most  of  which  have  yet  to  be  seen  in  ordered  intermetallic  systems  and  most 
of  which  are  vibrationally  unstable.  The  only  phase  which  we  found  that 
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Table  3:  Structural,  energy,  and  vibrational  stability  results  for  the  structures  of  IrV 
discussed  in  this  paper.  All  calculations  use  the  QE  pseudopotentials  discussed  in  the 
text.  “Source”  indicates  the  origin  of  the  unit  cell,  either  from  experiment  or  an  unstable 
phonon  in  the  indicated  structure.  “Atoms”  is  the  number  of  atoms  in  the  primitive  cell. 
Lattice  constants  are  given  in  the  standard  crystallographic  convention,  e.g.  the  lattice 
constants  of  the  full  orthorhombic  cell  are  given  for  the  Cnunm  and  Cmcrn  space  groups.  All 
primitive  cells  have  a  =  7  =  90°.  The  rows  “Ir”  and  “V”  give  the  Wyckoff  positions  of 
the  atoms.  “Volume”  is  the  minimum  energy  volume  of  the  structure  per  formula  unit,  in 
(A3).  “Energy”  is  the  energy  of  the  structure  below  the  Llo  structure,  in  meV/formula 
unit.  _ 


Structure 

Llo 

Cmnmi 

Cmcm 

Source 

Exp. 

U3(0H)L10 

phonon 

a3(00|) 

Cmmm  phonon 

Space 

Group 

P4  /  mmm 

Cmnmi 

Cmcm 

Atoms 

2 

8 

32 

a  (A) 

2.762 

7.260 

7.198 

b  (A) 

2.762 

11.110 

11.166 

c(A) 

3.668 

2.774 

11.136 

/? 

90 

90 

90 

Ir 

(la)  (000) 

(4e)(H0) 

(4g)(.267  0  0) 

(4d)(ii0) 

(4e)(.269  0  0) 
(4g)(.273  .254  ±) 
(4g)(.268  .002  4) 

V 

(Id)  (HO) 

(4j)(0  .364  i) 
(4j)(0  .115  1) 

(4f)(0  .358  .134) 
(4f)(0  .632  .117) 
(4f)(0  .123  .132) 
(4f)(0  .895  .117) 

Volume 

27.97 

27.97 

27.97 

Energy 

0 

10.98 

20.61 

Stability 

Unstable 

Unstable 

Unknown 
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Table  4:  Continuation  of  Table  [3l 


Structure 

GaaPts 

Amm2 

Imm2 

Source 

Ss(i^)  Ll0 

phonon 

Pj  Cmmm 
phonon 

Hi(i00) 

Amm2  phonon 

Space 

Group 

Cmmm 

Amm2 

Imm2 

Atoms 

8 

8 

32 

a  (A) 

7.222 

2.788 

7.176 

b  (A) 

7.800 

7.146 

11.144 

C  (A) 

3.974 

11.230 

11.190 

p 

90 

90 

90 

Ir 

(4e)(H0) 

(4h)(.273  0  i) 

(4d)(0  .230  .752) 
(4d)(0  .270  .998) 

(4c) (.256  0  .252) 
(8e)(.268  2.50  .252) 
(4c)(0.223  0  .752) 
(4c) (.256  0  .998) 
(8e)(.268  .250  .998) 
(4c)(.223  0  .498) 

V 

(2a)  (000) 

(2b)(|i0) 

(4j)  (0  .218  I) 

(2b) (|  0  .352) 

(2b) (i  0  .625) 

(2b)  (i  0  .125) 

(2b) (i  0  .899) 

(4cl)  (0  .131  .352) 
(4d)(0  .381  .355) 
(4d)(0  .132  .625) 
(4d)(0  .384  .625) 
(4d)(0  .130  .125) 
(4d)(0  .380  .125) 
(4d)(0  .131  .898) 
(4d)(0  .381  .895) 

Volume 

27.98 

27.97 

27.96 

Energy 

20.25 

24.15 

25.12 

Stability 

Unstable 

Unstable 

Unknown 
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Table  5:  Continuation  of  Tables  IHlmd  HI 


Structure 

Cm 

Rh5V3 

adrV 

Source 

Ai(OiO) 

Amm2 

phonon 

r4  GasPts 
phonon 

Exp. 

Space 

Group 

Cm 

Amm2 

Cmnnn 

Atoms 

32 

8 

4 

a  (A) 

13.303 

4.070 

5.816 

b  (A) 

11.147 

7.029 

6.767 

c(A) 

6.652 

7.822 

2.823 

P 

65.3 

90 

90 

Ir 

(2a) (.999  0  .504) 

(2a) (.265  0  .974) 

(2a) (.485  0  .534) 

(2a) (.755  0  .994) 

(2a) (.870  0  .256) 

(2a) (.140  0  .715) 

(2a) (.360  0  .276) 

(2a) (.626  0  .745) 

(4b) (.992  .249  .521) 
(4b) (.260  .252  .983) 
(4b) (.365  .748  .267) 
(4b)(.133  .251  .729) 

(4d)(0  .219  .724) 
(4e)(|  .281  .001) 

(4j)(0  .216  1) 

(4b)(.178  .130  .344) 

V 

(4b) (.175  .380  .362) 
(4b) (.312  .134  .625) 
(4b) (.312  .385  .625) 
(4b) (.064  .125  .117) 
(4b) (.060  .375  .133) 
(4b) (.447  .130  .906) 
(4b) (.450  .380  .888) 

(2a)  (0  0  .032) 

(2a)  (0  0  .457) 

(2b)  (i  0  .218) 

(2b) (i  0  .793) 

(4g)(.297  0  0) 

Volume 

27.97 

27.97 

27.78 

Energy 

25.67 

41.14 

67.25 

Stability 

Unknown 

Metastable 

Ground  State 
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is  metastable  is  the  experimentally  observed  RI15V3  phase,  albeit  in  a  50-50 
composition. 

However,  none  of  the  unstable  structures  shows  any  sign  of  relaxing 
into  the  ground  state  structure.  In  fact,  looking  at  Tables  |3]|5l  we  see 
that  all  of  the  structures  except  adrV  have  approximately  the  same  volume, 
27.97  A3/formula  unit  in  IrV.  The  ground  state  cHrV  structure,  one  the  other 
hand,  has  a  somewhat  smaller  volume,  27.78  A3/formula  unit.  While  this 
may  seem  insignificant,  the  volume  change  from  one  phase  to  another  is  a 
signal  of  a  first-order  transition,  although  only  weakly  first  order  here. 

This  paper  also  shows  a  mechanism  for  generating  new  candidate  in- 
termetallic  phases:  look  for  vibrationally  unstable  modes  in  a  high  energy 
structure,  construct  a  supercell  which  will  mimic  that  mode  within  the  frozen- 
phonon  calculation,  and  relax  the  cell.  Using  this  method  we  found  five  new 
structures,  as  well  as  one,  RI15V3,  which  had  been  seen  before  but  which  we 
had  not  been  aware  of  until  this  research  started.  Many  more  new  structures 
can  undoubtedly  be  derived  from  just  this  system,  but  further  research  in  this 
area  is  currently  restricted  by  the  time  needed  to  search  the  Brillouin  zone 
of  a  given  crystal.  For  the  eight  atom  supercells  computing  the  phonon  spec¬ 
trum  at  45  q-points  with  a  reasonable  number  (<  64)  of  processors  took  days 
to  weeks  to  complete,  depending  on  the  symmetry  of  the  crystal.  Calculating 
the  phonons  for  larger  unit  cells  will  of  course  take  longer,  although  we  will 
have  fewer  q-points  to  consider,  alleviating  some  of  the  0[iV3]  increase  due  to 
the  larger  unit  cell.  One  of  onr  research  goals  will  be  to  construct  a  “set-and- 
forget”  mechanism  which  will  search  the  Brillouin  zone  of  an  initial  structure 
(here  Ll0)  find  all  unstable  modes  with  supercells  containing  a  given  number 
of  atoms,  relax  those  modes,  and  repeat,  until  all  such  structures  have  been 
found.  Even  restricting  ourselves  to  all  possible  binary  intermetallicsjl]  this 
will  require  a  Grand  Challenge  computational  program. 
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